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ABSTRACT 

Soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) are thought to be magnetars: 
neutron stars powered by extreme magnetic fields. These rare objects are characterized by repeated 
and sometimes spectacular gamma-ray bursts. The burst mechanism might involve crustal fractures 
and excitation of non-radial modes which would emit gravitational waves (GWs). We present the 
results of a search for GW bursts from six galactic magnetars that is sensitive to neutron star f- 
modes, thought to be the most efficient GW emitting oscillatory modes in compact stars. One of 
them, SGR 0501+4516, is likely <~ 1 kpc from Earth, an order of magnitude closer than magnetars 
targeted in previous GW searches. A second, AXP IE 1547.0—5408, gave a burst with an estimated 
isotropic energy > 10 44 erg which is comparable to the giant flares. We find no evidence of GWs 
associated with a sample of 1279 electromagnetic triggers from six magnetars occurring between 
November 2006 and June 2009, in GW data from the LIGO, Virgo, and GEO600 detectors. Our 
lowest model-dependent GW emission energy upper limits for band- and time-limited white noise 
bursts in the detector sensitive band, and for /-mode ringdowns (at 1090 Hz), are 3.0 x 10 44 e^ erg 
and 1.4 x I0 47 dl erg respectively, where d\ = and d 50i is the distance to SGR 0501+4516. 

These limits on GW emission from /-modes are an order of magnitude lower than any previous, and 
approach the range of electromagnetic energies seen in SGR giant flares for the first time. 



Search for Gravitational Wave Bursts from Six Magnetars 



3 



Subject headings: gravitational waves — stars: magnetars 



1. INTRODUCTION 



Magnetars are isolated neutron st ars (NS) powered by 
extreme m agnetic fields (~ 10 15 G) (Duncan & Thomp- 
The magnetar model explains the observed 



son 



1992) 



properties of two classes of rare objects, the soft gamma 
repeaters (SGRs) and the anomalous X-ray pulsars 
(AXPs): compact X-ray sources with long rotation peri- 
ods and rapid spindowns which sporadically emit short 
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0.1s) bursts of soft gamma rays (for a review see 



Mereghetti||2008| ) Fewer than twenty SGRs and AXPs 
are known. The total isotropic burst energies rarely ex- 
ceed 10 42 erg. However, three extraordinary "giant flares" 
(GFs) have been observed in ^30 years from SGRs in 
our galaxy and the Large Magellanic Cloud: one from 
SGR 0526—66 in 1979 with a n observed total isot ropic 
energy of ~ 1.2 x 10 44 d§ 5 erg (Mazets et al. 1979); one 
from SGR 1900+14 in 1998 with 4.3 x lO 44 ^ erg (jTanaka 



10 46 d 2 5 erg 



Terasawa et al. 


20051 


so evidence that 


some 



in 2004 with - 5 x 
where d n = d/(nkpc). There is a 
short gamma ray bursts (GRBs) were in fact extragalac- 
tic GFs. GRB 070201 might have been a GF located 
in the Andro meda galaxy with an isotropic energy o f 



1.5 x 10 45 erg flMazets et al. 



2008 



Abbott et al. 



2008a); 



and GRB 051103 might haveTieeh a GF in M SI with an 
energy of 7.5 x 10 46 erg ( |Frederiks et"al]|2007| ). 

Although still poorly understood, magnetars are 
promising candidates for the first direct gravitational 
wave (GW) detection for several reasons. First, a sud- 
den localized energy release could excite non-radial pul- 
sational NS modes. Bursts may be caused by untwisting 
of the global interior magn etic field and associated crack - 
ing of the solid NS crust ( Thompson & Duncan 1995), 
or global reconfiguration of the internal magnetic field 
and a ssociated deformation of the NS hydrostatic cquilib 



rium ( |Ioka|2001| |Corsi k. Owen"|2011| ). The lowest- order 
GW emitting mode, the /-mode, is damped principally 
via GW emission and would ring down with a predicted 
damping time of 100-400 ms and with a frequency in the 
1-3 kHz range dependin g on the nuclear eq uation of state 
and NS composition (IBenhar et al. 20041) , putting these 



signals in the band of interferometric GW detectors (see 
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Figure!!]) . Second, precise sky locations and trigger times 
from electromagnetic (EM) bursts allow us to reduce the 
false-alarm rate and increase sensitivity relative t o all- 



sky all-time searches such asAbadie et al. (20101. Fi- 
nally, magnetars are among the closest of potential GW 
burst sources. 

GW signals from magnetars would give us a new win- 
dow through which to probe the stellar physics and struc- 
ture. However, quantitative predictions or constraints 
on the amplitude of GW emission associated with mag- 
netar bursts are relatively few and highly uncertain (see 



eg.JIoka|2001[ |Owen l 2005l |Horowitz fc Kadau|2 009 ; Corsi 
fc Owen||2011[ |Kasniyama fc: loka||2011| |Levin &: van 
Hoven|2011 "Jjnence it is not clear when we might begin 
to expect a detection. 

It may turn out that the magnetar burst mechanism 
does not excite global NS /-modes. If the outburst dy- 
namics are confined to surface layer modes, the crust 
torsional oscilla tions might emit GWs at frequencies of 
- 10 - 2000 Hz ( |McDermott et al.|1988r i. It is also possi- 
ble that although the crust is a plausible site f or trigger 
ing, b ursts are confined to the magnetosphere(Lyutikov 



2006 1 , although even in this case /-modes might be ex- 
cited either directly or via crust/core hydromagnetic cou- 
pling. Finally we note that it is not yet clear if GFs and 
common bursts are caused by the same mechanism. The 
lack of theoretical understanding underlines the impor- 
tance of observational constraints on GW emission. 

We present results from a search for GW bursts associ- 
ated with magnetar EM bursts using data from the sec- 
ond ye ar of LIGO's fifth science run (S5y2, Abbo tt et al. 



2009a); Virgo's first science run (VSR1, |Accrncse et al. 
2008[); and the subsequent LIGO and GEO astrowatch 



period (A5), during which the principal goal was detec 
tor commissioning, not data collection. The S5y2 epoch 
involved the three LIGO detectors: a 4 km interferome- 
ter in Louisiana and two interferometers (4 km and 2 km) 
in Washington. The VSR1 epoch added the Virgo 3 km 
detector to the global network. The A5 epoch included 
onl y the LIGO 2 km detector and the GEO 600 m detec - 
tor ( |Grote fc the LIGO Scientific Collaboration! |2010[ ). 
The Virgo and GEO600 detectors are located in Italy 
and Germany, respectively. 

This is the third search for GWs from magnetars sen- 
sitive to /-mode ringdowns. The first ( |Abbott et al. 
2008b| ) included the 2004 SGR 1806-20 GF, a 2006 
storm of bursts from SGR 1900+14, and 188 other 
events from SGRs 1806—20 and 1900+14 occurring be- 
fore November 2006. Upper limits on /-mode GW en- 
ergy emission at 1090 Hz ranged from 2.4 x 10 48 erg to 
2.6 x 10 51 erg, and upper limits on band- and time- 
limited white noise bursts at 100-200 Hz ranged from 
3.1 x 1 45 erg to 7.3 x 10 47 erg. The second ( [Abbott et al 
2009b) focused on the 2006 SGR 1900+14 storm, "stack- 
ing" GW data correspo nding to individual b ursts in the 
storm's EM lightcurve (IKalmus et al.||2009|. An upper 



limit on /-mode emission at 1090 Hz of 1.2 x 10 erg per 
burst was set on a stack of the 11 brightest storm bursts, 
an order of magnitude lower than the unstacked limit on 
the storm. 

During the S5y2, VSR1 and A5 epochs of the search we 
present here, 1217 soft gamma-ray bursts from six mag- 
netars were listed by the interplanetary network of satel- 
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LIGO Hanford 4km 2007-03-18 
LIGO Livingston 4km 2007-08-30 
LIGO Hanford 2km 2007-05-14 
Virgo 3km 2007-09-05 
GEO 600m A5 2008 




frequency (Hz) 



Fig. 1. — Best detector noise spectra from the LIGO and Virgo 
detectors during S5/VSR1 and the GEO600 detector during A5. 
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Fig. 2. — Each mark represents a burst from one of the six 
magnctar sources. Exceptional events are annotated in the figure; 
SGR 0501+4516 and SGR 0481+5279 were discovered in the A5 
epoch. The LIGO S5yl epoc h was the subject of the first /- mode 
search ( Abbott ct al. 2008b|. The current search includes bursts 
which occurred during tne LIGO S5y2 and Virgo VSR1 epochs 
for which usable data were available, as well as the A5 astrowatch 
commissioning period. The VSR1 epoch, which is a subset of the 
S5y2 epoch, is indicated by cross-hatching and darker shading. 
(NB: unabbreviated source names are given in Table[TJ) 



lites ^Jj or IPN (TablefTjand Figure[2j). Four of the sources 
are being examined for GW signals for the first time. 
Two of those (SGR 0501+4516 and SGR 0418+5729) 
are thought to be much closer to Earth than SGRs ex- 
amined in previous GW searches. SGR 0501+4 516 might 
be associated with the supernova remnant HB9 ( |Gaensler 
fc Chatter jcc 2008), which is (800 ± 400) pc from Earth 
( jLeahy e^l 'ian 2007); proper motion measurements could 



exclude this association. The probable locations of both 
SGR 0501+4516 and SGR 0418+5729 in the Perseus 
arm of our galax y imply distances of ~1— 2 kpc ( |van der 
Horst et aL|'2Q10] ). AXP IE 1547.0-5408 (also known as 
SGR 1550-5418) gave two exceptional bursts on 2009 
January 22. Observations of expanding rings around 
the source, caused by X-ray scattering off dust sheets, 
set the source distance at 4-5 kpc and imply an EM 
energy for one or both of these "ring event" bursts of 
10 44~45 erg (|Tiengo et al.||2010|, comparable to the GFs. 



In addition to the 1PJN triggers, we include eight triggers 
from the Fermi GBM detector: seven bright AXP IE 
1547.0-5408 bursts and one SGR 0418+5729 burst. We 
also identified 54 individual peaks in a storm from SGR 
1627-41 lasting ~ 2000 s by combining the 15-2 5 keV 



<)<X !)!) 



and 25-50 keV Swift/BAT 64ms-binned light curve; 
and selecting peaks above 450 counts / 64 ms. The search, 
thus includes a grand total of 1279 EM triggers. 

2. METHOD 

W e analyze magnetar bursts using the strategy 
from Abbott et al. (2008b), which is less dependent on 
a part icular emissio n m odel than the stacking approach 
of Abbott et al. ( 2009b[). The analysis i s performed b y 
the Flare pipelme fKalmus et aLl|2007| |Kalmus|[2008] ), 
which produces a time-frequency excess power pixel map 
from calibrated detector data streams in the Fourier ba- 
sis. Pixels are characterized by excess power relative to 
the background ("loudness") and loud adjacent pixels are 
grouped into "events." The generalized pipeline accepts 
arbitrary networks of GW detectors by including detector 
noise floor measur ements and an tenna responses in the 
detection statistic (| Kalmus | |2010[ ). We divide the search 
into three frequency bands: 1-3 kHz where /-modes are 
predicted to ring; and 100-200 Hz and 100-1000 Hz. We 
include the latter two frequency bands in order to search 
also at lower frequencies where the detectors are most 
sensitive (see Figure]!]). Although there are no pre- 
dictions of GW burst signals from magnetars at these 
lower frequencies, we note that quasiperiodic oscillations 
(QPOs), lasting for tens of seconds and possibly associ- 
ated with stellar torsional modes, have been observed in 
GF EM tai ls at frequencies as low as 18 Hz and as high 
as 1800 Hz flStrohmayer fc Watts | 2005| jlsrael et al.||2005 ' 



Steiner fc Wa tts 2009). QPOs in the tai l of the aJ0~ 
GF w ere targeted by a tailored GW search) Abbott et al 
2007) distinct from the one presented here. 

Xs in Abbott et al. ( 2008b ) , we choose 4 s signal regions 
centered on each EM trigger time. Delays bet ween EM 



http://ssl.berkeley.edu/ipn3 



and GW em ission are unlikely to be significant (Kalmus 
et al. 2009[ ); the 4 s duration accounts for uncertainties 
in the geocentric EM peak time due e.g. to satellite trig- 
gering algorithms and rounding. Overlapping signal re- 
gions are merged. We analyze 1000 s of background on 
either side of each signal region (2000 s total) in order 
to estimate the significance of events in that signal re- 
gion. Background regions are not necessarily continuous, 
as we require the same detector network coverage and 
data quality as for the signal region; in addition, signal 

http://heasarc.gsfc.nasa.gov/FTP/swift/data/ 
obs/2008_05/00312582000 

http: //heasarc. gsfc.nasa.gov/FTP /swift / data/obs / 
2008_05/00090056009 
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TABLE 1 

Summary of source sky locations and estimated distances. The nominal distances d N 



FOR SETTING UPPER LIMITS. ENERGY UPPER LIMITS CAN BE SCALED TO ANY DISTANCE d VIA THE FACTOR d 2 / d\, . SOME EM TRIGGERS 



OCCURRED WHEN THERE WAS NO GW DATA AVAILABLE (i.E. N = 0). 



ARE THE DISTANCE USED IN THE SEARCH 
JV' 



Distances (kpc) EM triggers 

Estimated Nominal total 



Source 



SGR 0418+572^] 



Position 

J2000 



Analyzed with JV detectors 

N = 1 N = 2 N > 3 



SGR 0501+451^] 



04 h 18 m 33.867± 0.35 s 
+57°32'22.91 ± 0.35" 



05 h 01 m 06.8± 1.4 s 
+45°16'35.4± 1.4" 



-2, 0.8±0.4 



1 



166 



105 



21 



-540^1 15 h 50 m 54.11± 0.01 s 
-54°18'23.7± 0.1" 
16 h 35 m 51.84±0.2 s 
-47°35'23.31±0.2" 



AXP IE 1547.0 



SGR 1627-410 



4-5, 9, 4 4 

n±o3 it 



844 
56 



315 512 



-1.5' 



SGR 1806-2! 



SGR 1900+l|] 



18 h 08 m 39.32±0.3 s 
-20°24'39.5±0.3" 



8.7 



6.4-9.8 



10 



207 



11 



:-!<> 



136 



19 h 07 m 14.33± 0.15 s 
+09°19'20.1±0.15" 



3-9, 12-15 



10 



a position: Woods ct al. (2009); distance: Esposito ct al. (2010); van dcr Horst ct al. (2010) 

b position: [Evans fe O sborne ( 2 008 [ ); dist ance: [van der Jiorst|et al. (2010); Ga ensler fc Cnatte rjcc ( 20081 ; [Leahy fc Tian 
c position: |Oarniio et al.|lpiUUY[l: distance: |Tiengo"et ai.| j^uiujl ; [Cainiio et'al.l | [^uu/[ |; |Geltand fc <jaensierJ ( |^uuY| i 
d position: Wachtcr ct al. (2004); distance: Gorbcl et al. (1999) 



(2007) 



e position: 
f position: 



rvaplan et al.| [I 



IVail et al. (1999) 



CT20U2[I ; distance: bibby ct"aT 1 20TI5 F]" Cameron et al.| (|2005[ ) 
1995 fTdistance: Marsden efaf |2D0Tl; Vrba et al. ] |2o00" 



regions of other magnetar bursts are masked out. Sig- 
nal and background regions are chosen after data quality- 
cuts have been applied to the GW data, so as to remove 
data segments coincident with instrumental or data ac- 
quisition problems, or excessive noise due to challenging 
environmental conditions. For the S5y2 portion of this 
search, we applied category 1 and 2 data quality cuts 
(i.e. cutting only data certain to be unfit for analysis) as 
described inAbadie et al. (2010). For A5, which focused 
on detector commissioning, the boolean "science mode" 
designator and other basic data quality treatments were 
applied to the data, but the full categorical data quality 
treatment was not performed. Statistically significant 
events in the signal regions from any epoch are subject 
to follow-up investigations before being considered de- 
tection candidates. Follow-ups might include correlation 
with environmental data channels and more refined esti- 
mates of significance. 

We set model-dependent upper limits on /-mode ring- 
downs with circular and linear polarizations and frequen- 
cies sampling the range for /-modes (1-3 kHz, which ac- 
counts for plausible NS equations of states and magnetic 
fields), and with a decay time constant of r — 200 ms. We 
observed no more than 15% degradation in strain upper 
limits using ringdowns with t in the range 100-300 ms 
as compared to the nominal value of 200 ms. We set ad- 
ditional limits on band- and time- limited white noise 
bursts with 11ms and 100 ms durations (motivated by 
observed rise times and durations of magnetar burst light 
curves) spanning the 100-200 Hz and 100-1000 Hz search 
bands. While these frequencies are chosen principally to 
explore the detectors' most sensitive region below the /- 
mode frequencies, the observed range of QPO frequencies 
provides astrophysical motivation. Upper limits depend 
on the frequency sensitivity of the detectors (Figure [l}. 

Simulations are constructed using knowledge of the tar- 
get ma gnetar's sky location and the EM burst time. Fol- 
lowing 



h 2 



Abbott et al. 

J — oo +, 



(2008b), h 



S2 

TSS 



h 2 

"TSS+ 



dt and h + , x (t) are the two GW po- 
larizations. The relationship between the GW polariza- 



tions and the detector response h(t) to GW signals ar- 
riving from an altitude and azimuth (0, (f) and with po- 
larization angle if) is: 

h{t) = F+{e,<f>,i>)h+{t) + f x (M,# x W, (i) 

where F + (9, cf>, ip) and F x (8, <fi, -0) are the antenna func- 
tions for the source at (9,<j>). The polarization angle for 
each simulation was randomly chosen from a flat distri- 
bution between and 2ir. The GW emission energy (if 
the integrand is averaged over inclination angle) is 



E, 



GW 



= 4ird 2 



16nG 



dt. 



(2) 



We estimate model-dependent upper limits on -Egw or 
h rss for a given signal region as follows: 

(1) We determine the loudest event in the signal region. 

(2) For a specific simulated signal type, we inject a sim- 
ulation at a specific £qw & n d h rss in a randomly 
selected 4 s interval of the background data and find 
events in that region. We compare the loudest sig- 
nal region event to the loudest event with a cluster 
centroid time near the known injection time (within 
100 ms for ringdowns and within 50 ms for white noise 
bursts) . 

(3) We repeat ^ for a range of £"gw & n d ft-rss values, and 
at each value we determine the fraction of injections 
with associated events louder than the loudest signal 
region event. 

(4) We repeat ^ using different simulated signal types. 
For each signal type, we estimate the 90% detection 
efficiency loudest event upper limit, E^Sff or h^°, at 
which 90% of injection events would be louder than 
the loudest signal region event. 

3. RESULTS AND DISCUSSION 

We find no evidence of a GW signal in any of the sig- 
nal regions analyzed. The loudest event of the search 
occurred at 2009 January 22 05:48:43.2 UTC and was 
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the only event with a false-alarm rate below our prede- 
termined follow-up threshold of 1/(3 x 4808 s) = 6.9 x 
1CP 5 Hz as estimated via extrapolation from the 2000 s 
local background region. This event cannot be consid- 
ered a GW candidate because it was found when only 
the Hanford 2 km detector was observing and was co- 
incident with a strong glitch caused by fluctuations in 
the AC power picked up by a magnetometer, and thus is 
highly likely to be an instrumental artifact. 

We estimate /i^ s % and E^S^ for each signal region, 
which depend on detector sensitivities and antenna fac- 
tors, the loudest signal region event, and the simulation 
waveform type. -Eqw u PP er limits also depend on nom- 
inal source distance and can be scaled to any source 
distance d via the factor (d/d-^) 2 . 

Figure[3] shows -Eqw upper limits for each of the EM 
triggers from the six magnetar candidates, for each wave- 
form type. The complete table of upper limits is on- 
linJ™| We spotlight bright bursts from SGR 0501+4516 
and AXP IE 1547.0-5408; however it is unknown 
whether E EM and E GW are correlated. Table[2] presents 
Eqw and h ISS upper limits for three exceptional EM trig- 
gers; and for a burst from SGR 0501+4516 occurring in 
the signal region centered at 2008 August 23 16:31:22 
UTC which yielded the lowest limits of the search. Each 
was analyzed with a network of the LIGO 2 km and 
GEO600 detectors. The SGR 0501+4516 burst with the 
larges t EM fluence (2.21 x 10~ 5 erg/cm 2 ( Aptekar et al 



(1) Additional GFs could push down upper limits on 7 



p90% / 771 



2009 1 , which corresponds to a 1 kpc isotropic energy of 
2.7 x 10 39 erg) occurred in the signal region centered at 
2008 August 24 01:17:58 UTC. The two candidate pro- 
genitor bursts for the expanding X-ray rings around AXP 
IE 1547.0-5408 occurred at 2009 January 22 6:45:14 
UTC and 6: 48:04 UTC. with e stimated isotropic E EM of 
10 44-45 erg( |Tiengo et al.|2()l()[). Table[2]also gives upper 
limits on the ratio 7 = Eq ^ /Eem for the three bursts 
with £em estimates. The 7 upper limits for the two ring 
bursts were estimated using Eem = 10 45 erg, and beat 
the best prev ious upper limits on 7, set for the SGR 
1806-20 GF ( | Abbott et al.|2008b| ), by a factor of a few. 

Superscripts in Table[5] give uncertainties at 90% con- 
fidence. The first is uncertainty in detector amplitude 
calibrations. The second is the statistical uncertainty 
(via the bootstrap method) from using a finite number 
of injected simulations. Both are added linearly to final 
h ISS upper limit estimates; corresponding uncertainties 
are added to -Egw upper limit estimates. 

Our best -Egw /-mode upper limits are an order of 
magnitude lower (better) than the best /-mode lim- 
its from previous searches, and approach the range of 
EM energies seen in SGR GFs for the first time. The 
best SGR 0501+4516 /-mode limit of 1.4 x 10 47 erg (at 
1090 Hz and a nominal distance of 1 kpc) probes be- 
low the available energy pre dicted in a f racti on of the 
parameter space explored in Ioka (2001) and |Corsi & 



Owen (2011), the predicted maximum being ~ 
10 4a erg. The best 100-200 Hz white noise burst limit 
of 3.5 x 10 44 erg is — for the first time — comparable to 
the -Eem seen in "normal" GFs. 

Improved upper limits and perhaps detection will come 
in the future via the following routes: 

100 https:/ /dcc.ligo.org/cgi-bin/DocDB/ShowDocument?docid= 



(2) An analysis which stacks isolated bursts (from e.g. 
SGR 0501+4 516 and AXP IE 154 7.0-5408) using 
the method o fj Abbott et"al] ( |2009b[ ). Stacking 100 or 
more bursts observed with a constant detector sen- 



sitivity, as in Abbott et al. (2009b), might yield up 
to an additional order of magnitude improvement in 

7790% 

(3) The GW detectors will become more sensitive. Sec- 
ond generation detectors (Advanced LIGO and Ad- 
vanced Virgo) are expected to begin observing by 
2015, promising more than two orders of mag- 
nitude improvement in Eqw sensitivity over the 
LIGO 2 km + GEO600 network w hich obser ved SGR 
0501+4516 JAbbott et al.||2009a[ ). Recently [LevTiTfc 
van Hoven ( 201 ip made semi-quantitative predictions 
on /-mode excitations in GFs. Their predictions are 
pessimistic that an /-mode signal from a GF at 1 kpc 
would be detectable in the second generation, though 
they do not consider crustal cracking. Third gener- 
ation detectors could yield two additional orders of 
magnitude in energy sensitivity. 

We look forward to further predictions on GW emission 
amplitudes from these enigmatic sources. 
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Fig. 3. — ^gw u PP er limits for the entire SGR burst sample for various circularly/linearly polarized ringdowns (RDC/RDL) and white 
noise burst (WNB) signals (see Section)^. For each of twelve waveform types, we show six rows of dots marking upper limits for the sources 
(from top to bottom) SGR 1900+14 (violet); SGR 0418+5729 (purple); SGR 1627-41 (orange); SGR 1806-20 (green); SGR 0501+4516 
(teal); and AXP IE 1547.0-4508 (blue) for that waveform type. The limits shown in Table[2]for SGR 0501+4516 and AXP IE 1547.0-4508 
are indicated in the figure by circles. 
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TABLE 2 

gw strain and energy upper limit estimates at 90% confidence (/l 9 ,° s % and e 9 -, ^), for the burst trigger yielding the lowest 
■z^gw upper limits (top left), the brightest sgr 0501+4516 burst (top right), and the two 'ring' events from axp ie 
1547.0—5408 (bottom). Upper limits on the ratio 7 = Bq^/Eem are given when estimates for Eem are available; for the 
ring events 7 = -Eq^/IO 45 erg. Upper limits were estimated using the circularly and linearly polarized ringdowns 
(RDC/RDL) and white noise burst (WNB) waveforms (see Section]!}. Uncertainties, from detector calibration and using a 
finite number of injected simulations, are added to the final upper limit estimates. These are given for the h Isa limits as 
superscripts, with the first showing detector calibration uncertainty and the second showing statistical uncertainty 

from finite injected simulations. 
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